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Abstract 

We propose a method to determine the quantum numbers, which we call the rigged 
configurations, for the solutions to the Bethe ansatz equations for the spin-1/2 isotropic 
Heisenberg model under the periodic boundary condition. Our method is based on the 
observation that the sums of Bethe’s quantum numbers within each string behave 
particularly nicely. We confirm our procedure for all solutions for length 12 chain 
(totally 923 solutions). 
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1 Introduction 


Bethe’s seminal solution to the isotropic Heisenberg model under the periodic boundary 
condition [B] in 1931 is one of the prototypical theories on the quantum integrable systems. 
Basic procedure of the algebraic Bethe ansatz is as follows (see, e.g., BEEll). First, we 
solve the set of algebraic equations called the Bethe ansatz equations (see equation flT^ in 
the main text). Next, by using the solutions to the Bethe ansatz equations, we construct 
eigenstates of the Hamiltonian (see equation ffTT]) b The main problem which we will consider 
in the present paper is to show that the whole procedure involved in the Bethe ansatz method 
is mathematically well-dehned. We pursue the problem in the same setting as was treated 
by Bethe. 

For a long time it has been observed that there are several subtle points about the 
procedure. One of such obstructions is the problem called the singular solutions (see equation 
firsji b Recently we have much progress on this issue (see |AVl IBMSZl INWl IHNSl IKSldal 
IKS 14b] 1 and now we attain fairly good understanding of the problem. However, the very 
structure of the solutions to the Bethe ansatz equations itself still remains in rather foggy 
situation. 

In the analysis of the solutions to the Bethe ansatz equations, it has been customary to 
assume that roots of the solutions take a particular form called the strings. This idea was 
already apparent in the original Bethe’s paper. However, as we will explain shortly after, the 
notion “string” has rather elusive nature and for a long time it has been a difficult problem 
to understand how to utilize the string structure in a proper manner. One of the well-known 
attempts on this problem is Takahashi’s theory pp. Let us explain in more detail. Let N be 
the length of the spin chain and let £ be the number of down spins. Then a main assumption 
in the derivation of Takahashi’s quantum number is to suppose that the strings take exactly 
the following formji] 


a + bi, a + {b — 1)*, a + {b — 2)i, ..., a — bi, (a G M, 6 G Z>o/2). (1) 


However, when £ > 3, the assumption ([T]) has serious difficulty. Indeed, as we can easily 
see in examples, the real part a as well as the intervals between successive roots are not 
unique if the lengths of the strings are larger than 2. As the result, Takahashi’s quantum 
numbers are not uniquely defined and also they are not half-integers. In a previous paper 
|S15] . we proposed to seek an alternative to Takahashi’s quantum numbers. In particular. 


in that paper we showed that the correct quantum number for the exceptional real solution 
[EKSj is different from the one derived by Takahashi’s quantum numbers. 

In the present paper, as a continuation of a previous work |S15] . we propose a method 
to assign quantum numbers to strings of roots. For this purpose, we start from the so-called 
Bethe’s quantum numbers (see Section[3]). Here, roughly speaking, Bethe’s quantum numbers 
appear as phase factors after taking logarithm of the Bethe ansatz equations. By dehnition, 
Bethe’s quantum numbers are uniquely dehned and exactly integers or half-integers. Then 
our basic observation is that the sum of Bethe’s quantum numbers associated with all roots 


^We want to point out that the assumption of the validity of the string shape 0 allows to guess some 
explicit expressions for q-weight multiplicities (Kostka-Foulkas polynomials) appearing in the representation 
theory of Lie algebras of type A. All these formulas have been proven rigorously in [KRj and lead to 
applications in combinatorics, representation theory and discrete integrable systems. 
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of a given string behaves in a particularly simple manner. This is a remarkable property 
since individual Bethe’s quantum numbers behave in a rather complicated way. Based on 
this observation we propose a method to determine complete set of the quantum numbers, 
which we call the rigged configurations, from Bethe’s quantum numbers. We confirmed these 
observations for all solutions of iV = 12 case by using the numerical data given in |GD1 IHNS] . 

As we noted previously, we do not have thorough understanding of the string pattern 
which appear in the solutions to the Bethe ansatz equations. However, our main aim in the 
present paper is to propose a method which could be made mathematically rigorous once we 
obtain sound understanding of the string structure. Our expectation relies on the fact that 
Bethe’s numbers are uniquely defined and exactly half-integers. 

To our best knowledge, Bethe’s quantum numbers had been introduced and studied in the 
original paper by Bethe [B], formulas (37a) and (37b), and has been studied more recently 
in ra in their study of the set of real solutions to the Bethe ansatz equations for i = 2, 
and in [HCj concerning the deviation of solutions to the Bethe ansatz equations from having 
exact string structure with numerical evidence for N = 10. 

The outline of the present paper is as follows. In Section [2l we provide necessary facts 
about the Bethe ansatz method and the rigged conhgurations. In Section [3] we collect 
necessary facts about Bethe’s quantum numbers. In SectionlUwe provide the main algorithm. 
We provide examples in Section[5l Finally we conclude in Section ^ Besides the main article, 
there are supplementary tables (see remarks after Conjecture [3]) o 


2 Algebraic Bethe ansatz 

In this section we briefly overview some basic facts concerning the spin-1/2 isotropic Heisen¬ 
berg model, also known as the XXX model. The space of states of the spin-1/2 XXX chain 
of length N is the complex space , which one can identify with the tensor product of N 
copies of C^, namely 


N 

i=i 


( 2 ) 


The Hamiltonian describing the interaction of particles has the form 

7 ^ 

«" = 7 - I")' (3) 

^ k=l 

where we assume = a", a G {x, y, z}, and 


C7 


X 





(4) 


^ After completing the present work, we noticed that the paper |DG15) appeared. In this paper the 
authors also study Bethe’s quantum numbers. However our work is significantly different from theirs since 
our main motivation is to understand the situation which appear when £ > 3 whereas their paper considers 
the case £ = 2 exclusively. Note that by a result of m the real part “a” of the string 0 is uniquely defined 
when £ = 2 (except for the exceptional real solutions). 
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stand for the Pauli matrices, and 




a 

k 


Let us introduce the vectors 


v+ = 


and V- = 


and the global vacuum (ground state) as 

|0)jv = u+ (8) 

Let us consider the local transfer matrix 


v+ e s^N- 


LkW = 


AIjv + 


2^k . 


i(^k 


where ± ia^, and consider the monodromy matrix 


(5) 

( 6 ) 

(7) 

( 8 ) 


Ta,(A) = L«(A)Ljv_i(A)...Li(A)= ( S(A))^ (”) 

Recall that the local transfer matrices satisfy the relation R 12 L 1 L 2 = L 2 L 1 R 12 for some 
matrix R G 0 where R satishes the quantum Yang-Baxter equation 


R 12 R 13 R 23 — -R23-Ri3-Ri2- (10) 

The quantum Yang-Baxter equation flTOl) implies that the transfer matrices tr |c 2 T/v(A) and 
tr|c 2 T/v(p) (Cq is the auxiliary space) commute for different values of A and /r. As a con¬ 
sequence of the relation (1T0|) . we can deduce that [R 7 v(A),-Bv(h)] = 0 for different values of 
the parameters A and /i. The key observation due to Bethe, reformulated in the language of 
the algebraic Bethe ansatz, is that the vector 


\k7v(Ai, . . . , Xe) — Bn{Xi) ■ ■ ■ i?7v(A£)|0)Ar 


( 11 ) 


is an eigenstates of the Hamiltonian if and only if the parameters Ai,..., A^ satisfy the 
following system of algebraic equations, known as the Bethe ansatz equations (BAE(£) for 
short) 


Afc + 
Afc — 


N 


j-j- Afc - Xj 


+ i 


J = 1 




(t = i,■■■,<). 


( 12 ) 


We are interested in the so-called physical solutions to BAE(£). The set of physical 
solutions is divided into the regular solutions, that is, the solutions (Ai,..., A^) which do not 
contain a pair (Aq, A^) = (|, —|), and the set of physical singular solutions of the form 


A = 



(13) 
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In the case of the regular solutions, one has 



j=l j 4 


However for singular solutions, the energy is not defined. It was suggested to “re¬ 
solve” singularity of by using a deformed form of singular solutions, namely. 



( 16 ) 


and consider the limit 



( 16 ) 


It was conjecture by |NW] . that a 7 ^ 0 if and only if 



(17) 


and 



j=3 2 


We call the singular solutions which satisfy the condition f[T7)) the physical singular solution. 


The corresponding energy is |KS14b] 



It is also conjectured [HNS] that the set of physical singular and regular solutions to BAE(£) 


exhausts the set of all physical solutions to BAE(£). 

Finally, let us define the rigged configurations. The rigged configurations (z/, I) are com¬ 
prised of a partition u = {ui,... ,ui), called the configuration, and a sequence of non-negative 
integers Ij, called the riggings, satisfying the following condition. Define the vacancy number 
Pk{i^) by the following formula 


Pk{i^) = N - 2(z/j H- \- ul) = N - 2 ^ min(fc, Uj) 


(19) 


where Uj represent components of the transposed partition u'. If we regard the rigged con¬ 
figuration as the set of pairs {v, I) = {(z^i, Ji),..., (z/;, Ii)}, then the riggings must satisfy 


0<lj< Pu^{iy) 


( 20 ) 
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for all j = 1,..., / . We note that in the rigged configuration theory, the order of the riggings 
associated with rows of the conhguration v with the same length is not essential. Therefore 
if we have the subset {(r'fc, • • •, Ik,m)} within (i/, J), we assume that Ik,i < • ■ ■ < Ik,m 

to erase ambiguity. 

In our earlier paper |KS14a] we conjectured that there is a one to one correspondence 
between the set of physical solutions to BAE(.^) and the set of the rigged configurations (z/, I) 
where the Young diagrams representing the partitions u have i cells. Let us define the flip 
operator on the set of physical solutions to BAE(£) as follows 

K : {Ai, A 2 ,..., Af} I —> {—Ai, —A 2 ,..., —Af}. (21) 

Then another conjecture in |KS14aj states that the corresponding rigged configuration are 
connected by the following operation on the rigged configuration (z^,/). Namely, we replace 
each element (z/j, Ij) by (z/j, — Ij) and reorder the riggings if necessary. This property 

is useful in the following discussion. 

3 Bethe’s quantum numbers 

In what follows, we use function Arctan(z), defined as an analytic continuation of the function 
arctan(x) {x G M, arctan(O) = 0) to all complex plane with branch cuts (i, +ioo) and 
{—i,—ioo) along with the imaginary axis. For example, Arctan(a:i) = iarctanh(a:) and 
Arctan(—xi) = —i arctanh(a:) for x> 1. If we go across the branch cuts, we have 

Arctan(xi) — hmArctan(xi — e) = tt, 

Arctan(—xi) — lim Arctan(—xi + e) = —vr 

for X > 1 and e > 0. 

We use the following formula 

log —j-T = Arctan(z) + {2n + l)7ri (n G Z). (22) 

We verify the formula as follows. We see that the differentiations of the both sides give 
In order to determine the remaining constant, we compare the values of the both sides at 
z = 0 as follows; 

log = log(— 1 ) = log(e’^* • = (2n + l)7ri (n G Z), 

0 — -z 

2i Arctan(O) = 0. 
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3.1 Regular case 

First let us consider a regular solution. By using the formula fl2^ . we take the logarithm of 
the Bethe ansatz equations flT^ : 


log 


Afc + 


1 N N I 

2 


n 

j=i 


Aj % 


e 

= —N ■ 2i Arctan(2Afc) + 2i Arctan(Afc 


\j) + (—A^ + £ — l)7ii + 2mri (n e Z). 


By the Bethe ansatz equations, this should coincide with log(l) = 2m7ri for m G Z. Writing 
2Jk = 2n — 2m — A^ + £ — 1, we obtain 


Jk 


N 


/ 

2 Arctan(2Afc) 

V 


^ ^ Arctan(Afc 

3¥=k 



(23) 


We call the half integers Jk the Bethe’s quantum numbers. 

Here we describe the basic property of the Bethe’s quantum numbers in relation with the 
multiplication by ( — 1). 


Proposition 1. Suppose that we have the following two solutions of the Bethe ansatz equa¬ 
tions 


{Ai, A2,..., Xi}, 

{Ai, A2,..., A^} := {—Ai, —A2,..., —Xe}. 

Let Jq, (resp. Ja) he the Bethe’s quantum number corresponding to Aq (resp. X^). Then we 
have Ja = — Ja- 

Proof. Since Arctan(— 2 ;) = — Arctan( 2 ;), we obtain the result by fl23l) . ■ 

3.2 Singular solutions case 

Now let us consider the singular solution flT^ when N is even. The computation for Jk 
{k > 3) is same with the previous case. For the computations of A and J 2 , we should be 
reminded that the function Arctan(z) has the logarithmic singularities a.t z = ±i. Recall 
that we assume e G M in ffTSj) . We introduce the deformation of singular solutions 

HO) = -l + 
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where r G M>o and study the limits of 


N 

2 Arctan(2Ai(0)) — — Arctan(Ai( 6 ') — ^2(9)) 

N 

2 Arctan( 2 A 2 ( 6 *)) — — Arctan(A 2 ( 6 *) — Ai( 6 *)) 


when r —)■ 0. As a result we obtain the Bethe numbers Ji{0) and J 2 {G) which depend on 
the choice of 9. The following two graphs show Ji{0) (left) and J2{9) (right) for the solution 
{i/ 2 , -i/ 2 } oiN = 12 . 



We note that the sum Ji{9) + J2{9) is constant for —| < 9 < ^. 

4 Main algorithm 

We start from a given physical solution A = (Ai,..., A^) to the Bethe ansatz equations. We 
divide A into subsets called strings. Here we say that S = {rji,... ,rii}, {j]j = aj + ibj) forms 
a string if 

ttj - ttj+i Ri 0, bj - bj+i Ri 1 (j = 1,1). 

By a result of Vladimirov one can assume that = rj; if / > 1. We define the real part 
of the string by Re(S') = Re^rji). 

We take a component of the solution A to the Bethe ansatz equations with maximal 
imaginary part and starting from that component we form a string of maximal length. 
To proceed, we erase the string created, and apply the above procedure for the remaining 
components of the solution A. As a result we partition all component of the solution A into 
a collection of distinct strings. 

One of the main ingredients of our algorithm is the sum of the Bethe numbers corre¬ 
sponding to a string S 

J{S) = (24) 

rjGS 

A striking feature of the quantity J{S) is that it exhibits very stable behavior even if indi¬ 
vidual Bethe numbers behave in a complicated way. 

Example 2. Let us think about solutions of V = 12 which contain a 1-string and a 5-string. 
We arrange the solutions according to the real parts of the 5-strings. We depict the solutions 
on the complex plane with the Bethe’s quantum numbers for each root. 












#1 , #2 , #3 



• 7/2 

. 13/2 

. 13/2 


, 11/2 

. 

, 9/2 

- 5/2 , 

_ 5/2 - 3/2 . 

, 1/2 . -1/2 , 

, - 3/2 


• 13/2 

• 9/2 

. 11/2 


• 9/2 

. 15/2 

, 15/2 

#4 

#5 

#6 



13/2 

9/2 - 15/2 , 

t 


9/2 

11/2 - 11/2 

1 

1/2 , 

- 3/2 3/2 

- 3/2 3/2 

. - 1/2 


11/2 

- 11/2 - 9/2 , 



15/2 

- 9/2 - 13/2 ^ 


#r 

#8 

#9 


- 15/2 

- 15/2 , 

- 9/2 , 


- 11/2 

- 9/2 

-13/2 , 


3/2 

. 1/2 . - 1/2 . 

. 3/2 ^ -5/2 , 

. 5/2 

- 9/2 

- 7/2 , 

- 11/2 , 


- 13/2 

- 13/2 , 

-7/2 . 



Let S be the 5-string of each solution. Then for solutions #1,... ,#4, we have Re(S') > 0 
and J{S) = 45/2. For solutions # 6 ,... ,#9, we have Re(S') < 0 and J{S) = —45/2. Note 
that individual Bethe’s quantum numbers behave in a rather complicated way. ■ 

In solution ^5 of the above example, we have the roots ±a (a G M, a 1) near the 
origin. Although we cannot decide which root should belong to the 5-string, this ambiguity 
does not affect the following procedure. See solutions 7^6 and #16 of Section E] for another 
type of such ambiguity which does not affect the procedure also. 

4.1 Configuration 

Starting from the solution A, we divide it into strings Si,Sp. Let the lengths of Sk be 
Ik- Then the partition z/ which appears in the rigged conhguration (z/, I) is u = (R,..., lp)~^. 
Here the symbol (/i,..., denotes the decreasing order of numbers (/i,..., #). 
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4.2 Riggings 

We determine the rigging corresponding to a string in the following way. If a string S satishes 
Re(S') = 0, then the corresponding rigging is the half of the corresponding vacancy number. 
Otherwise consider the set of all solutions with a conhguration u, denoted by 

Sol(z/) = = {A(“)}. 

To start with we compute all Bethe numbers for Sol(i^). 

The main tool to construct bijection between the set of solutions to the Bethe ansatz 
equations and rigged conhgurations is the following string crossing rule. Let us consider a 
string S satisfying Re(S') > 0 and length n. Let ..., Si^ be strings such that 

• length of = k < n, 

• Re(^) < Re(^iJ 

ioT k = 1,... ,q. Similarly let ..., Sj^, be strings such that 

• length of Sj^ = l'^.> n, 

. Re(SjJ < Re(S) 

for A; = 1,..., g'. Then we consider 

J(S) := J{S) - A(4,?7,) + ^ A(/fc,n), A{l,n) := I + n - 3. (25) 

k=l k=l 

Let us £x some z/* of the partition u and denote the multiplicity of Ui by m. Let be 
the strings of length Ui within A^"^ Here we assume that the strings of the same length Ui 
satisfy Re(^2) > Re(^2) > ''' > Dehne 

= max{J(^£5^)}. 

OL ’ 

Suppose that we have Re(S'|“^^) > 0. Then the corresponding rigging is 

P„[v) - {M„., - J(S5)}. 

If we have Re{S^°‘\) < 0, we can use the flip operation to obtain the rigging. 

Then our main conjecture is the following. 

Conjecture 3. The above procedure gives one-to-one correspondence between the set of 
physical solutions to the Bethe ansatz equations and the rigged configurations. 

We conhrm the conjecture for the length 12 chain based on the numerical data by EH. 
As supplementary data, we provide plots of the solutions to the Bethe ansatz equations for 
A = 12 case (except for simpler ones). The solutions are arranged according to the real part 
of the largest string. For each solution, we assign the numbers corresponding to the solution 
numbers in |GDj . 
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Remark 4. Roughly speaking, the riggings represent positions of strings. Namely the 
larger rigging corresponds to the larger real part Re(S'). In |S15j . we observed that in several 
examples we can readily read off the riggings if we arrange the solutions according to the 
real parts of the largest strings. We conhrm that this algorithm works for all solutions of 
length 12 if the number of the strings is at most 3. In other cases, the appearance of the 
solutions in the ordered set of solutions becomes more complicated. Nevertheless we observe 
there is a strong relation between the riggings and the positions of the strings even in such 
case. ■ 

4.3 Explicit formulas for A4 

We expect that the numbers M depend only on the length N of the chain and the shape of 
the conhguration z/. We make the following conjecture. If z/2 < 3, we have 

_ (AT - zyj -max{l,z/^})z/i 
1 ^ 1,1 2 

If z/i > 1 and z /2 = 1, we have 


+ z/i + z /2 — 2z/j — 3 

M„ =-^- 

We checked these formulas for length 12 chain. 


5 Examples 

We consider the case of z/ = (3, 2,1) and N = 12. Below we give a list of the solutions 
depicted on the complex plane which we arrange according to the real part of the 3-string. 
To each root we provide the corresponding Bethe numbers. Following the diagrams we give 
tables of the explicit numerical values which we need in the computation of J. Note that the 
solution 7^11 contains a singular string. See also |S15l Section 5] for an alternative method 
to determine the riggings. 


Group 1. 


#1 #2 


-7/2 

• 

7/2 

• 

.5/2 

-7/2 

• 

7/2 

• 

5/2 

-5/2. 


-3/2 


-5/2 

• 

-5/2 

• 


9/2 


9/2 


#3 #4 


-7/2 

• 

.7/2 

.5/2 

-7/2 

• 

.7/2 

5/2 

-.1/2 


• 

1/2 

-5/2 

• 

-5/2 

• 


9/2 


9/2 
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#5 #6 #7 


-7/2 

• 

7/2 

• 

5/2 

-7/2 

• 

‘ 9/2 

• 

,3/2 . 

-7/2 

• 

‘ 9/2 

• 

.3/2. , 

• 

3/2 

• 

3/2 

• 

5/2 

-5/2 

• 

-5/2 

• 

-5/2 

• 


9/2 


11/2 


11/2 


#1 

0.54241927 

-0.45810568 + 0.50017785* 

0.54455699 + 0.99639165* 
-0.45810568 - 0.50017785* 

0.54455699 - 0.99639165* 
-0.71532188 

#2 

0.52708058 

-0.64127830 + 0.50335013* 

0.52957875 + 0.99660493* 
-0.64127830 - 0.50335013* 

0.52957875 - 0.99660493* 
-0.30368149 

#3 

0.49893578 

-0.69284388 + 0.50515234* 

0.50200196 + 0.99724969i 
-0.69284388 - 0.50515234* 

0.50200196 - 0.99724969* 
-0.11725196 

#4 

0.46564665 

-0.71976299 + 0.50614240* 

0.46941736 + 0.99809739* 
-0.71976299 - 0.50614240* 

0.46941736 - 0.99809739* 
0.035044606 

#5 

0.42430010 

-0.73869344 + 0.50683156i 

0.42960641 + 0.99941330* 
-0.73869344 - 0.50683156* 

0.42960641 - 0.99941330* 
0.19387395 

#6 

0.38490522 + 0.01906127* 
-0.75221326 + 0.50729383* 

0.36730804 + 0.99179719* 
-0.75221326 - 0.50729383* 

0.36730804 - 0.99179719* 
0.38490522 - 0.01906127* 

#r 

0.23056669 

-0.76056174 + 0.50745313* 

0.23083274 + 0.99967059* 
-0.76056174 - 0.50745313* 

0.23083274 - 0.99967059* 

0.82889133 


Group 2. 


#8 

#9 

#10 

#11 



.9/2 

• 11/2 

► 11/2 

,9/2 


. 7/2 

,5/2 

,5/2 



. 1/2 . 

t 1/2 ^ 

1/2 . 

-1/2 . 

-5/2 

-3/2 

. -1/2 


1/2 


7/2 

5/2 

*5/2 



11/2 

* 13/2 

13/2 

' -9/2 


#12 #13 #14 
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#8 

0.20669577 

0.10578435 + 0.50000000* 

0.20597572 + 1.00038608* 
0.10578435 - 0.50000000* 

0.20597572 - 1.00038608* 

-0.83021590 

#9 

0.059726272 

0.10847310 + 0.50000000* 

0.06007063 + 0.99927337* 
0.10847310 - 0.50000000* 

0.06007063 - 0.99927337* 
-0.39681373 

#10 

0.010757119 

0.06941354 + 0.50000000* 

0.01249979 + 0.99958901* 
0.06941354 - 0.50000000* 

0.01249979 - 0.99958901* 
-0.17458378 

#11 

0.018539900* 

0.50000000* 

0.99377501* 

-0.50000000* 

-0.99377501* 

-0.018539900* 

#12 

-0.010757119 
-0.06941354 + 0.50000000* 

-0.01249979 + 0.99958901* 
-0.06941354 - 0.50000000* 

-0.01249979 - 0.99958901* 
0.17458378 

#13 

-0.059726272 
-0.10847310 + 0.50000000* 

-0.06007063 + 0.99927337* 
-0.10847310 - 0.50000000* 

-0.06007063 - 0.99927337* 
0.39681373 

#14 

-0.20669577 

-0.10578435 + 0.50000000* 

-0.20597572 + 1.00038608* 
-0.10578435 - 0.50000000* 

-0.20597572 - 1.00038608* 

0.83021590 


Group 3. 


#15 


#16 


#17 


#18 


-11/2 


-11/2 


-9/2 


-9/2 


• 

5/2 

• 

5/2 


5/2 


5/2 

-3/2 


-3/2, 




-5/2. , 


-5/2 

• 

-3/2 

• 

-3/2 

• 

-1/2 

• 

• 

7/2 

• 

7/2 

• 

7/2 

• 

7/2 

-9/2 

#19 
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#15 

-0.23056669 

0.76056174 ± 0.50745313* 

-0.23083274 ± 0.99967059* -0.23083274 - 0.99967059* 
0.76056174 - 0.50745313* -0.82889133 

#16 

-0.38490522 ± 0.01906127* 
0.75221326 ± 0.50729383* 

-0.36730804 ± 0.99179719* -0.36730804 - 0.99179719* 
0.75221326 - 0.50729383* -0.38490522 - 0.01906127* 

#17 

-0.42430010 

0.73869344 ± 0.50683156* 

-0.42960641 ± 0.99941330* -0.42960641 - 0.99941330* 
0.73869344 - 0.50683156* -0.19387395 

#18 

-0.46564665 

0.71976299 ± 0.50614240* 

-0.46941736 ± 0.99809739* -0.46941736 - 0.99809739* 
0.71976299 - 0.50614240* -0.035044606 

#19 

-0.49893578 

0.69284388 ± 0.50515234* 

-0.50200196 ± 0.99724969* -0.50200196 - 0.99724969* 
0.69284388 - 0.50515234* 0.11725196 

#20 

-0.52708058 

0.64127830 ± 0.50335013* 

-0.52957875 ± 0.99660493* -0.52957875 - 0.99660493* 
0.64127830 - 0.50335013* 0.30368149 

#21 

-0.54241927 

0.45810568 ± 0.50017785* 

-0.54455699 ± 0.99639165* -0.54455699 - 0.99639165* 
0.45810568 - 0.50017785* 0.71532188 


Let us explain how our algorithm works for the examples listed. 

3-strings. For solutions #1, ..., ^5, we have J = 21/2 since there are no shorter strings 
on the right of the 3-string. 

In solution 7 ^ 6 , we have two roots 0.38 ± 0.02L We classify one of them to 1-string and 
the other to 3-string. In any case, 1-string is on the right of the 3-string. Thus we have 
J = 23/2 — A(3,1) = 21/2. Similarly we have J = 21/2 in solution ^7. 

In solution 7 ^ 8 , we have J = 21/2 since there are no shorter string on the right. 

In solutions ^^9 and ^^10 we have J = 25/2 — A(3, 2) = 21/2 since the 2-string is on the 
right of the 3-string. 

In solution 7 /=!!, the real part of the 3-string is 0. Thus the corresponding rigging is the 
half of the corresponding vacancy number P 3 ( 3 , 2 , 1 ) = 0 , that is, the rigging is 0 . 

In solutions #12, ..., #21, we have J < 0. 

Summarizing, we have AI 3 = 21/2. Since we have J = ±21/2 (except for the solution 
# 11 ), we conclude that the corresponding rigging is always 0 . 

2-strings. In solutions #1, ..., #7, #12, #13, and #14, the real part of the 2-strings are 
negative. 

In solution # 8 , we have J = 14/2 since the 3-string is on the right and the 1-string is on 
the left. 

In solutions #9 and #10, we have J = 10/2 ± A(2, 3) = 14/2 since the 3-string is on the 
left of the 2 -string. 

In solutions #15, ..., #21, we have J = 12/2 ± A(2, 3) = 16/2 since the 3-string is on 
the left of the 2-string. Note that the position of the 1-string does not affect J since we have 
A(2,l) = 0. 

To summarize, we have #^2 = 16/2. Thus the riggings for the two strings are 0 for 
solutions in Group 1, 1 for solutions in Group 2 and 2 for the solutions in Group 3. 
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1-strings. By a similar argument we see that Aii = 7/2. Then the riggings for each group 
is 1, 2,..., 6 from the smaller solutions number to the larger solutions number. 

The resulting rigged conhgurations agree with the ones obtained by a geometrical argu¬ 
ment |S15j or by the computation of Takahashi’s quantum numbers |GD] (see also [DG14j ). 

Remark 5. The above example reminds us about the box-ball system (see [SI2] for a review). 
The box-ball system is a discrete soliton system where solitons are “crystal” analogue of 
magnons. One of the important properties is that the rigged configurations provide action 
and angle variables of the box-ball system. In this picture each row z/j of the configuration v 
represents a soliton of the same length. If we look at the above examples, it is tempting to 
consider them as a propagation of the string of solutions where scattering of strings yields 
a phase shift A{m,n) = m + n — 2. Note that the phase shift in the case of the box-ball 
system is min(m, n) for the scattering of lengths m and n solitons. It will be interesting 
to remark that in the box-ball system the crystal analogue of the transfer matrices, which 
provides the quantum integrability of the box-ball system, generates the corresponding rigged 
configurations which are dynamical variables if we regard the box-ball system as classical 
integrable system |S07j . ■ 

6 Conclusion 

Summarizing, we describe an algorithm which associate collection of quantum numbers or 
riggings to solutions to the Bethe ansatz equations for the spin-1/2 Heisenberg model, start¬ 
ing from the set of Bethe’s quantum numbers corresponding to those solutions and its string 
structure for general N. We confirm our algorithm for all solutions to the Bethe ansatz 
equations for TV = 12. 
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